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Abstract. The unmatched X-ray resolution of Chandra allows probing the gas flow near qui- 
escent supermassive black holes (BHs). The radius of BH gravitational influence on gas, called 
the Bondi radius, is resolved in Sgr A* and NGC 3115. Shallow accretion flow density profiles 
n oc r"^ with /3 = 0.7 - 1.0 were found for Sgr A* and NGC 3115 with the help of Chandra. 
We construct self-consistent models with gas feeding and dynamics from near the Bondi radius 
to the event horizon to explain the observations. Gas is mainly supplied to the region by hot 
colliding stellar winds. Small-scale feedback such as conduction effectively flattens the density 
profile from steep /3 = 1.5 in a Bondi flow. We further constrain density and temperature pro- 
files using the observed radio/sub-mm radiation emitted near the event horizon. We discuss the 
present state of our numerical model and its qualitative features, such as the role of the galactic 
gravitational potential and the random motion of wind-emitting stars. 
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1. Introduction 

As indicated by both the theory and the observations a typical active galactic nucleus 
(AGN) is not particularly active ( Ho 2008) . The median Ed dington ratio in the Palomar 
survey of nearby AGNs is A = i^boi/^Edd ^ 10~^ (|HoII2009I ). Palomar surve y still misses 
AGNs like Sgr A*, which has an even lower ratio A ~ 10~^ ( Narayanl 12002). Extensive 



radio observations help to identify weaker AGNs. While Sgr A* is an estabhshed radio 
emitter, the activity in the n uclei of, e.g., M31 and NGC3115 is barely detected with 
the Jansky Very Large Array (ICrane et al.lll99 2: Wrobel & Nyland, 2012, subm.). AGNs 
have a short duty cycle (Gre ene fc H0II2OQ7I ) and are in a low-luminosity state most of the 



time. AGNs e nter a short bright quasar phase, when galaxy mergers feed the associated 
BHs with gas ^Hopkins et al.' 200 61 l2008h. AGNs shut off quickly after the mergers due to 



the large-scale feedback (Hopkins fc Hernquist 2009) . The material in an isolated galaxy 
can accrete in various ways onto the central BH, which may power Seyfert galaxies (e.g., 
Hopkins fc H ernQuistl l2006l ). When the gas s upply from other sources is small, the BH 



can be fed mostly by stellar winds (iHopki ns fc Her nquisti 12005 IHoII2009 I). Typical mass 



loss rates from nuclear star clusters lead to very sub-Eddington accretion rates. Colliding 
stellar winds produce a hot tenuous medium, which might have a long cooling timescale. 
In the a bsence of cooling the accretion occur s as a radiatively inefficient accretion flow 
(RIAF) (jNaravan et al.lll998l : lQuataertir200l[ ). Modeling BH feeding from stellar winds 
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and the associated accretion flow in low-luminosity AGNs (LLAGNs) is the purpose of 
the present work. We discuss the features of the accretion flow, which delivers only a 
small fraction of available matter to the BH event horizon. We construct a semi-analytic 
radial model and apply it to Sgr A* and NGC3115. 



2. Supply and Accretion of Stellar Winds 




Figure 1. The scheme of accretion in LLAGNs. Colliding stellar winds provide matter into the 
feeding region. Some of that matter outflows from the region, while a fraction accretes. The 
heat from the inner flow leads to supervirial outer temperature and stalled accretion. Image of 
Chandra satellite was produced by Chandra team. 

The matter is supplied near the radius of BH gravitational influence, the Bondi radius 
rs = GMbh/c^^ by stellar winds. The mass loss rate rhej in winds can be estimated 
based on properties of the nuclear stellar popula tion. The predicted rhej I'ate depends 
on the age of the stellar population (Jungwiert et al.l 200 ll), the ini tial mass function, 



metallicity, and the stellar evolution model ( Leitner fc Kravtsovll2011 ). Approximate cor- 



relations exi st of mass loss rate with V-band or B-band luminosity (Ciotti et ah 1991 



IPadovani fc Matteucci 1993 : h3[2009) for old stellar populations. Such relations are rela- 



tively easy to apply, since surface brightness in these bands was d etermined in many galac- 
tic n uclei as a by-product of searching for supermassive BHs (Kormendv fc Richstond 
Il995h . In the case of Sgr A* the mass loss rate was directl y computed from line proper- 
ties of individual wind-emitti ng stars (jMartins et al. 



As simulations show (e.g., Cuadra et al. "2008"), the gas may exist in a quasi-steady 
state over long periods of time near the Bondi radius. Accretion onto the BH and outflow 
balance the injection of stellar winds. The gas energetics is determined by the velocity 
of colliding stellar winds, the gravitational potentials of the BH and the galaxy, and 
the energy flux from the inner accretion flow. The wind-producing stars move in the 
gravitational potential at a speed Vgt^ while expelling winds at relative velocity v^j. As- 
suming no preferred direction of orbital angular momentum of stars, the specific energy 
of colliding winds is < v"^ v^^ > /2, where <> denotes the ensemble average. The 
gravitational potential of the galaxy is essential for setting the quasi-stea dy gas density. 



The galactic potential overtakes the BH potential at r ~ ^Otb for Sgr A* (|Schodel et al 
[2003 ) and at r - 2rB for NGC3115 (Shcherbakov et al. 2012, in prep). Gravity cre- 
ates a potential barrier for the outflow to escape from the feeding region. This potential 
barrier can be substantially underestimated, if the galactic potential is neglected. The 
percentage of the outfl owing material may gre atly vary. The accretion rate onto Sgr 



A* is ~ 10 ^M0yr ^ ( Shcherbakov et al. 20121 ), while the cluster wind mass loss rate 
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lO-^Moyr-i (|Mart,ins et a ll 120071) is much higher. Lastly, the energy balance can be 



modified by the deposition of energy diffused from the inner accretion fiow. Outward 
energy fiux can fiatten the density profil e and shut off the accretion. Also it can help 
drive an outfiow from the feeding region. iBlandford fc BegelmanI (|l999h state that "the 



binding energy of a gram of gas at a few Tg drives off 100 kg of gas from 10^ rg" , where 
rg is a gravitational radius. The reason for such a dramatic infiuence can be understood 
from pressure balance. The balance of the gas pressure and the gravitational force reads 

1 dp d(nkBT) _ GMrUp 
n dr ndr ' 

where n is particle density. The density profile is usually approximated as a power-law 
n (X . As the temperature increases slightly due to the deposited energy, the density 
slope P gets correspondingly shallower. A 25% temperature increase causes /3 to drop 
from 1.5 to 1.0, which leads to 300 times lower accretion rate M for r^/^g = 10^. 
The origin of the energy fiux may vary depending on gas density. Turbul ent diffu 



sion/convection is traditionally i nvoked to account for the energy fiux (N aravan et al 



2000l : [Ouataert fc Gruzinovll2000l ). The fiow is intrinsically turbulent and the large-scale 



eddies effectively transport the energy outwards. Heat conduction takes over for turbulent 
diffusion at low gs :as density, when th e mean free path of electrons becomes comparable 



to the fiow size ( Sharma et al.ll2008l ). Conductive heat fiux is significant in turbulent 



magnetized fiows, where the random magnetic field does not effectively trap the elec- 
trons (Narayan fc Medvedev 2001; Ruszkowski fc Oh 2011). As the problem is complex, 
and encompasses a huge range of radii tb/t^ ^ 10^ — 10^ preventing the direct numer- 
ical simulations, we developed a semi-analytic radial model to quantify the accretion in 
LLAGNs. 



3. Inflow-outflow Model with Conduction and Applications 

The description of our first model can be found in IShcherbakov fc Bagano3 ( 2010l ). 
We compute the radial profiles of gas density p, electron temperature Te, and proton 
temperature Tp. We compile the radial profiles of mass and energy injection in winds. In 
the latest version of the model we include the random velocity of stars to compute the 
effective wind velocity. The gravitational potential is similarly improved to include the 
galactic potential. The electron temperature Tg deviates down from Tp in the inner fiow 
at about lO^rg, where the electrons become relativistic. On the way to the BH different 
species are heated by pdV work and by viscous energy dissipation. The fiow had zero 
angular momentum in the fir st model, which may be a good assumption for the large 
values of effective viscosity a phakura fc Sunvae^ I ll973[ ) or outside of the circularization 
radius. The transport of angular momentum and the associated heating are modeled 
self-consistently in the latest version. The energy is transported outwards by conduction. 
Conductive heat fiux is proportional to the gradient of Tg, while conductivity at each 
radius r is a small fraction of a product of that radius r by the electron velocity Ve- 

The hot accretion fiow onto a BH has a specific radiation signature, which was observed 
i n nearby LLAG Ns. The gas temperature in the feeding region is about T ~ 1 keV 
(|Pel legrini '2005) . as produced by the colliding winds. The hot gas radiates in X-rays 
via a variety of physical proces ses such as heavy ion line emission. We use the emission 
model based on ATOMDB v2.0 (|Foster et al.ll2012[ ) to compute the spectrum in collisional 
equilibrium for an assumed metallicity Z. We employ custom-generated response matrices 
and photoelectric absorption to simulate the spectra as seen by Chandra and to compare 
with the observed ones. The X-ray spectrum is dominated by heavy ion emission for 
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gas temperatures in 0.3 — 1 keV range, which makes the estimate of gas density heavily 
dependent on metahicity. Inner density and temperature can be determined based on 
radiation from hot inner flow. The plasma emits mainly synchrotron radiation at the 
electron temperature Tg ^ 10^^ K. Assuming that the radio emission at the highest 
observed frequency of radio/sub- mm bump comes from near the event horizon, we can 
estimate the electron temperature and density based on that frequency and correspondent 
flux. The slope between the inner and the outer flows is P = 0.80 — 0.90 for Sgr A* 
(jShcherbakov et al.ll2012h and /3 = 0.7 - 0.8 for NGC3115 (Shcherbakov et al. 2012, in 



prep), while /3 = 1 was found in the outer flow ( Wong et al.H2011 ) of NGC3115. 
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